UNCLASSIFIED 


AD  NUMBER 

AD 90 1620 

NEW  LIMITATION  CHANGE 
TO 

Approved  for  public  release,  distribution 
unlimited 


FROM 

Distribution  limited  to  U.S.  Gov't, 
agencies  only;  Test  and  Evaluation;  Jun 
72.  Other  requests  for  this  document  must 
be  referred  to  Director,  Naval  Research 
Lab.,  Washington,  D.  C.  20390. 


AUTHORITY 

NRL  ltr  dtd  9  Nov  1997 


THIS  PAGE  IS  UNCLASSIFIED 


mwm 


1221.1 


DR.  GURSKY 


PLEASE  REVIEW  NRL  REPORT  .7404  TO  SEE  IF  THE 


DISTRIBUTION  STATEMENT  CAN  BE  CHANGED  TO  PUBLIC 


'Vive 


CHUCK  ROGERS 


MARY  TEMPLEMAN  (202)767-3425 


Prepare  2  copies  of  this  route  sheet  and 
forward  ALL  copies  together  with  neces¬ 
sary  correspondence  and  other  documents. 


DATE  OF  MATERIAL 

09/11/97 


OHIG.  IDENT.  SYMBOL  {Mail  Room  Fill  In) 


5227 /MT 


IDATE  MAILED 


•PURPOSES 


1.  FOR  INFORMATION 

2.  FOR  APPROVAL 

3.  PREPARE  REPLY 

4.  PREPARE  ENDORSEMENT 

5.  FOR  NECESSARY  ACTION 

6.  FOR  SIGNATURE 

7.  RETAIN  ENCLOSURES 

8.  RETAIN  COPY 


FORWARDING  OF: 

NRL  REPORT  7404  FOR  REVIEW  OF  DISTRIBUTION 
STATEMENT . 


R/S  NO.  (Mail  Room  Fill  In) 


RQ-NRL  521S/6  .'Rev.  2-82) 


AD90  1620 


NRL  Report  7404 


The  Equations  of  Motion  of  a  Satellite  in 
a  Local  Vertical  Coordinate  System 


£vFreijlrk  k  W.  Raymond 

Satellite  7  eehniques  Hriuu  li 
Space  System*  Division 


June  1.  |972 


D  D  C 

rn|gF'~"'~  r~- 

in]  21  !97 Z 

uikiSEir 

B 


NAVAL  RESEARCH  LABORATORY 
Washington,  D.C. 


Distribution  limited  to  U  S.  tinvernment  Agencies  only .  test  and  evaluation.  June  1072.  cither  requests  for  this  document 
must  be  referred  to  the  Director.  Naval  Research  Laboratory,  'Washington.  D.C.  20 JIM 


■ywr 


",'VT1 


CONTENTS 


Abstract .  ii 

Problem  Status .  ii 

Authorization .  ii 

INTRODUCTION .  1 

THE  LOCAL  VERTICAL  COORDINA  TE  SYSTEM .  1 

THE  EQUATIONS  OF  MOTION .  1 

MULTIPLY  LAUNCHED  SATELLITES .  3 

ORBITAL  POSITIONING,  PARKING,  AND  STATION 
KEEPING  .  5 

REFERENCES .  9 


1 


ABSTRACT 


The  equations  of  motion  of  a  satellite  have  been  derived  in  a 
local  vertical  coordinate  system.  The  method  of  derivation  and  the 
form  of  the  equations  provide  a  simplified  means  of  visualizing  and 
planning  the  maneuvering  within  multiple  satellite  constellations. 
The  equations  are  also  useful  in  station  keeping,  rendezvous,  and 
collision  avoidance  considerations. 

PROBLEM  STATUS 

This  is  an  interim  report;  work  is  continuing. 

AUTHORIZATION 

NR  L  Problem  401-17 


Manuscript  submitted  February  22,  1972. 


THE  EQUATIONS  OF  MOTION  OF  A  SATELLITE  IN  A 
LOCAL  VERTICAL  COORDINATE  SYSTEM 


INTRODUCTION 

The  equations  of  motion  derived  in  this  report  have  been  previously  published  (1,2). 
However,  the  method  of  derivation  and  the  form  of  the  equations  presented  herein  should 
permit  a  more  intuitive  comprehension  of  the  relative  motions  of  orbiting  masses.  This  in¬ 
formation  should  be  especially  useful  to  |>crsons  concerned  with  the  problems  of  satellite 
rendezvous,  positioning,  and  collision  avoidance. 

The  equations  are  derived  and  presented  in  a  “local  vertical”  coordinate  system,  which 
is  defined  as  a  coordinate  system  rotating  at  an  angular  velocity  and  phase  very  nearly 
equal  to  that  of  the  orbiting  objects  under  consideration.  It  is  felt  that  this  coordinate 
system  provides  an  intuitive  reference  system  for  considering  the  relative  motion  of  orbit¬ 
ing  objects.  The  validity  of  the  equations  will,  of  course,  require  that  the  objects  con¬ 
sidered  remain  in  very  similar  orbits.  Thus,  the  usefulness  of  the  equations  is  restricted  to 
satellites  injected  into  nearly  identical  orbits  and  to  the  terminal  phase  of  rendezvous 
maneuvers. 


THE  LOCAL  VERTICAL  COORDINATE  SYSTEM 

The  local  vertical  coordinate  system  is  defined  such  that  the  z  a\i>  is  always  in  the 
direction  of  the  radius  vector  R  from  the  center  of  rotation  to  a  point  in  the  cluster  of 
orbiting  objects.  This  radius  vector  is  constrained  to  rotate  with  a  constant  angular  veloc¬ 
ity  gj  so  as  to  descrilte  a  Keplerian  circular  orbit;  thus,  the  square  of  the  angular  veloc¬ 
ity  must  bo  inversely  proportional  to  the  cube  of  the  radius.  The  direction  of  the  v  axis 
is  parallel  to  the  angular  velocity  gj,  and  the  direction  of  the  .v  axis  forms  a  right-hand 
coordinate  system.  (It  is  parallel  to  the  velocity  of  the  tip  of  the  radius  vector  R.)  The 
origin  is  chosen  to  be  the  tip  of  the  radius  vector  as  shown  in  Fig.  1. 


THE  EQUATIONS  OF  MOTION 

The  effective  force  on  a  point  of  mass  m  in  the  rotating  local  vertical  coordinate  sys¬ 
tem  described  above  is  the  sum  of  the  external  forces,  i.e.,  the  gravity  force,  the  centrif¬ 
ugal  force,  and  the  Coriolis  force.  Thus, 

F  =  mr  =  F(ext)  -  -7——  -  mu  X  [coX  (R  +  r)J  -  2m[uX  (R  +  r)],  (1) 

IK  +  rr 
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satellite  position  where  v  is  the  gravity  constant.  By  expanding  the  vector 
^  triple  product  and  |R  +  r|-3,  it  is  easily  shown  that 


r  =  a(ext)  +  (R  +  r)|td|2 


[*-(*♦  i 


2  2£\-3/2 

+  IR 1/ 


(R  *r) 


-  (oj  •  r)o>-  2(cj  X  r)  , 


where  a(ext)  =  the  acceleration  due  to  external  forces. 
Note  that  R  -  0  and  e/|/f  |3  =  |co|2. 

'Ilie  absolute  value  of  R  is  at  least  3500  miles;  thus, 
for  jr |  less  than  350  miles  it  is  certainly  reasonable  to 
neglect  second-order  terms  in  |r|/|R|.  With  this  approxi¬ 
mation,  the  equations  of  motion  linearized  in  terms  of 
the  x,  y,  and  z  com|K>nents  of  r  become 


\  • 

CENTER, 

.t  =  a,  -  2coi 

(3) 

OF  ) 

ROTATION  / 

/  7 

y  =  a  o  -  co-y 

(4) 

i  =  n;,  +  3co 2z  +  2co.x  , 

(5) 

Fir.  1— Tin*  rotating  local  vertical  coor¬ 
dinate  system  when*  a  >,  and  are  res|H*etively  the  .v,  y,  and  z 

components  of  acceleration  from  external  forces  and 
co  is  the  magnitude  of  the  orbital  angular  velocity  of  the  reference  coordinate  system. 

These  equations  lend  themselves  nicely  to  solution  by  l-iplace  transforms.  The  solu¬ 
tions  jure 


*  =  *o  - 


rif-  /2a-,  \  (2'z0  4a,\  , 


2a; ,  4  x  o\ 

— -  6zn  _  -  sin  cut 

co-  co  / 


Vi,  cos  c 0/  +  (1  -  cos  col)  +  sin  c 0? 

co-  CO 


-  1  *  z 0  cos  col  +  — 77  +  4cy  +  -  ( 1 

CO  V  CO-  CO  / 


—  cos  col ) 


2a  1  -’o 

•»  T 

co-  CO 

/ 


sin  c of 


-  (2io  -  ^ 


Sa\t  —  I -  +  3.vq  + 

\  cj 


-  ”  6c02q  “  4.VyjfOSCOl 


sin  col 
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y 
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sin  a >(  +  yo  cos  cot 


(w 

+  +  3w^o  +  2ioj  sin  cot  +  +  *o  j  cos  to?. 


(10) 

(ID 


where  xq,  yo,  and  z q  are  the  initial  coordinates  and  xo.  yo,  and  io  are  the  initial  veloci¬ 
ties.  Figures  2  and  3  illustrate  the  trajectories  for  various  initial  conditions.  Motion  in  the 
y  direction  is  not  illustrated  because  it  is  always  simple  harmonic  motion  which  is  easily 
visualized  as  superimposed  on  the  x-z  plane  motion. 


MULTIPLY  LAUNCHED  SATELLITES 


One  of  the  problems  associated  with  multiple  payload  launches  is  the  possibility  of 
collision  among  the  payloads.  It  is  generally  recognized  that  the  probability  of  collision  of 
separately  launched  satellites  is  extremely  small;  however,  in  multiply  launched  satellites, 


Z 


KiH.  £ -Elliptical  motion  of  a  satellite  having  only  initial  velocity 


Fig.  3— Cycloidal  motion  of  a  satellite  having  only  .v  initial  velocity 
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which  arc  placed  into  essentially  the  same  orbit  by  a  single  launch  vehicle  the  probability 
of  collision  becomes  significant.  By  judicial  selection  of  the  separation  velocities  of  the 
payloads  relative  to  the  launch  vehicle,  measures  can  be  taken  which  will  minimize  the  pos¬ 
sibility  of  collision.  One  obvious  measure  is  to  separate  the  payloads  so  that  drag  forces 
will  cause  them  to  continually  separate.  For  instance,  consider  two  satellites,  one  having 
a  higher  drag  force.  The  satellite  with  the  high  drag  will  lose  more  energy  and  spiral  into 
higher  velocity  orbits.  Thus,  to  avoid  the  possibility  of  the  higher-drag  satellite  overtaking 
and  passing  the  low-drag  satellite,  it  should  be  initially  separated  into  a  higher  velocity 
orbit. 

Another  consideration  is  that  it  is  possible  to  separate  satellites  from  the  launch  ve¬ 
hicle  into  orbits  which  have  identical  jicriods  and  differ  only  in  inclination  and/or  eccen¬ 
tricity.  This  could  result  in  periodic  collisions  at  the  orbital  frequency  and  obviously 
should  be  avoided.  This  problem  can  bo  analyzed  using  the  previously  derived  equations. 

Let  the  local  vertical  coordinate  system  have  its  origin  at  the  center  of  mass  of  the 
system  of  satellites  and  launch  vehicle.  In  this  coordinate  system  the  trajectories  of  satel¬ 
lites  impulsively  ejected  with  initial  velocities  a'o,  >’o>  ai'd  !lIV  easily  obtained  from  these 
equations. 

Consider  the  case  of  an  initial  velocity  purely  in  the  x  direction,  or  flight  line,  where 
>  0  =  io  “  0.  Then 


x  =  -3jco*  +  4xo  sin  tot  (12i 

z  =  -^-(1  -  cos  cot)  (13) 

to 


and 


y  =  0  . 

rrhe  resulting  trajectory  in  the  x  =  z  plane  is  a  cycloid  as  shown  in  Fig.  3. 
Eliminating  the  parameter  tot  between  the  .v  and  z  equations  results  in 

(,v  +  3.v0t)2  + 

16*o~ 

,  •> 


(U) 


(15) 


which  is  an  ellipse  centered  at  .v  =  -3.vo(  and  z  -  2.Vq/co  with  major  and  minor  axes 
a  =  4-vo/vo  and  b  -  2xo/ui. 

Thus,  the  cyloidal  motion  can  also  be  thought  of  as  an  ellipse  falling  behind  at  the 
rate  of  3.xof.  A  negative  initial  x  velocity  would,  of  course,  ca  se  the  ellipse  to  fall  ahead. 
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Consider  the  case  of  an  initial  velocity  only  in  the  z  direction: 


(1  -  cos  c of) 
CO 


(16) 


z 0  . 

z  =  —  sin  u)t 

OJ 


(17) 


y  -  o  . 


(18) 


Since  x  and  z  are  truly  parametric  in  col,  time  can  be  eliminated,  which  results  in  a 
closed  orbit  described  by 


2isJ*  *  Oi  .  ,  . 


4  ?  2 
C-J2 


il 

CO2 


(19) 


This  is  an  ellipse  centered  at  z  =  0  and  x  “  2io/u)  with  u  *  2i<)/u;  and  b  =  zq  'ui  as 
illustrated  in  Fig.  2. 

Thus,  if  perturbations  such  as  drag  are  neglected,  a  satellite  separated  with  only  the 
?  component  of  velocity  will  return  to  its  starting  point  once  per  orbit  and  collide  with 
the  launch  vehicle  and  any  other  satellite  separated  simultaneously  with  only  z  velocity. 

It  is  easily  seen  that  satellites  with  only  y  initial  velocity  oscillate  with  a  simple  harmonic 
motion  along  a  line  in  the  y  direction  with  tin*  orbiL.1  period. 

It  is  also  easily  deduced  that  a  combined  initial  velocity  of  y  and  z  (no  .v)  serves  only 
to  rotate  the  ellipse  previously  described  out  of  the  .v-r  plane.  Thus,  any  separation  of  sat¬ 
ellites  with  no  a:  component  of  initial  velocity  will  result  in  a  collision  with  the  launch  ve¬ 
hicle  and,  for  simultaneous  separations,  with  each  other.  Simultaneous  separations  having 
the  same  initial  x  components  of  velocity  will  result  in  collisions  at  the  orbital  period,  re¬ 
gardless  of  their  initial  y  and  z  components. 

Thus,  to  minimize  the  possibility  of  collisions,  the  drag  effects  must  Ik*  considered 
and  the  separation  velocities  chosen  so  that  an  .v  component  is  always  greater  than  zero, 
unequal  to  the  x  component  of  any  other  satellite.  If  there  is  a  zero  ,v  component,  it  must 
lx*  assured  that  the  launch  vehicle  is  pushed  out  of  the  way  with  an  additional  impulse. 

Figure  4  shows  the  trajectories  relative  to  the  center  of  mass  of  a  typical  multiple  pay- 
load  separation. 


ORBITAL  POSITIONING,  PARKING,  AND  STATION  KEEPING 

If  Laplace  transforms  tire  used  again  to  solve  the  equations  of  motion,  it  can  be  shown 
that  a  constant  acceleration  commencing  at  t  =  0  and  ending  at  t  =  T  will  cause  the  change 
in  motion  described  by 
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Fij;.  4  Trajrclorit's  nl  typical  multiply  lauiteltcil 
uilvllilcs  relative  tu  the  center  of  nutv* 
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A.V  -  -^ll!  -  3oi  TAf  -  (cos  cu(T  +  Af)  -  coscoAf]  (20) 

2  co  co  * 

2(i  o 

+  — —  (sin  to(T  +  At)  —  sin  coA t ) 
co^ 

Ay  **  -—£■  (cos  to(T  +  A f )  -  coscoAf]  (21) 

A z  s  — r.1-  ~  -A  [cos  co(7' +  Af )  -  coscoAf]  (22) 

coJ  to* 


2oi 

co2 


(sin  to(T  +  A/)  -  sin  coAf] 


A.v  =  -3a  iT  + 


1 sin  co(7’  ♦  A 0  ~  sin  coA<) 
co 


(23) 


+ 


— ^  (cos  co(7' +  Af)  “  eosuAf] 
co 


Av  =  —  (sin  oo(7'  +  AO  “  sin  coAf) 
co 


(24) 


=  (sin  o o(T  +  Af)  -  sin 

u) 


ooAf]  +  — ■  (cos  to (7  Af)  ■  cos  coAf  ]  .  (2o) 


whore  Af  is  measured  after  the  termination  of  the  tlmist. 


It  should  be  noted  that  if  the  thrust  period  T  is  an  integral  multiple  of  the  orbital 
period  2n/oj,  the  oscillating  terms  go  to  zero.  This  means  that  the  thrusting  did  not  change 
the  eccentricity  or  inclination  of  the  orbit. 


Usually  station  keeping  is  accomplished  with  thrusters  directed  along  die  flight  line 
(the  .v  direction).  The  purpose  is  to  control  A.v  and  A.v  since  these  coordinates  most  di 
rectly  affect  the  angular  position  of  the  satellite  in  orbit.  If  the  other  parameters  which 
are  usually  of  little  consequence  for  station -keeping  purposes  are  neglected,  the  station¬ 
keeping  equations  become 


A.v 


3tn7’2 

9 


3aj  7'Af 


A.v  =  —  3a  1 7'  . 


If  the  station-keeping  equations  are  divided  by  the  radius  of  the  orbit,  A x/li  and 
A x/li  become  an  angle  and  an  angular  rate  respectively,  thus  providing  a  good  approxima¬ 
tion  to  the  equations  necessary  for  the  spacing  and  parking  of  satellites  with  a  flight-line 
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thruster.  Satellites  can  also  be  positioned  with  r-axis  thrusting;  however,  the  rate  of 
change  of  x  exists  only  during  the  period  of  i  thrusting. 

The  Ax  for  each  orbital  period  of  i  thrust  is  4*03 /to2.  With  a  three-slug  satellite  and 
a  10-micropound  z  thrust,  A.v  per  orbital  period  of  thrust  is  40  ft.  A  10-micropuund  .v 
thrust  will  change  the  position  by  195  ft  during  one  revolution  of  thrusting  and  then  the 
Ax  will  continue  to  change  by  390  ft  each  revolution  thereafter. 

Although  z  axis  thrusting  is  not  very*  effective  and  cannot  in*  used  for  parking,  it  may 
be  useful  in  applications  requiring  very  fine  positioning.  It  does  have  the  advantage  of  not 
requiring  a  three-axis  stability. 

Considering  Kqs.  (20)  through  (25)  it  is  apparent  the  terminal  phase  of  a  rendezvous 
maneuver  consists  of  measuring  the  relative  positions  Ax,  Ay,  and  A z  and  relative  veloci¬ 
ties  Ax,  Ay,  and  A z  of  a  target  vehicle,  then  applying  thrust  so  as  to  make  the  relative 
position  and  velocity  zero. 


The  magnitudes  and  durations  of  the  thrusts  in  each  of  the  axes  provide  six  p.iram 
eters  that  can  be  controlled  to  optimize  a  rendezvous  maneuver.  Kquations  (20) through 
(25)  provide  the  basis  for  the  necessary  computations. 


Most  spacecraft  maneuvers,  other  than  station  keeping,  are  accomplished  with  short 
thrusts  from  relatively  large  motions.  They  can  be  treated  as  impulsive  velocity  changes. 
The  rendezvous  or  positioning  equations  can  thus  be  simplified  to 


Ax 


•  4xy  .  2c.> 

“d.Vi )t  f - Sill  uj/  -  -  <1  -  cuscOf) 

CO  CO 


A  vy  . 

Ay  =  - sin  cot 

CO 

A z  “  sin  cot  +  — ^  (1  -  eoscot) 
co  co 

Ax  ”  — 3xy  +  4 Xy  cos  cot  -  2c y  sin  cot 
Ay  *  yy  cos  cot 
A i  c  iy  cos  cot  +  2.vy  sin  c ot  , 

where  xy,  yy,  and  Cy  are  now  considered  impulsive  velocity  changes. 


(2l»  I 


(271 


(28) 


(29) 


(30) 

(31) 


These1  equations  have  proven  to  be  extremely  useful  in  visualizing  and  planning  ma¬ 
neuvers  necessary  to  position  satellites  in  orbit  and  provide  a  means  ot  predicting  the  com¬ 
plex  constellations  that  can  be  formed  with  multiple  satellites. 
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